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ABSTRACT 



Hadronic 7-ray emission from supernova remnants (SNRs) is an important tool to 
test shock acceleration of cosmic ray protons. Tycho is one of nearly a dozen Galactic 
SNRs which are suggested to emit hadronic 7-ray emission. Among them, however, it 
is the only one in which the hadronic emission is proposed to arise from the interac- 
tion with low-density (~ 0.3 cm -3 ) ambient medium. Here we present an alternative 
hadronic explanation with a modest conversion efficiency (of order 1%) for this young 
remnant. With such an efficiency, a normal electron-proton ratio (of order 10~ 2 ) is 
derived from the radio and X-ray synchrotron spectra and an average ambient den- 
sity that is at least one-order-of-magnitude higher is derived from the hadronic 7-ray 
flux. This result is consistent with the multi-band evidence of the presence of dense 
medium from the north to the east of the Tycho SNR. The SNR-cloud association, in 
combination with the HI absorption data, helps to constrain the so-far controversial 
distance to Tycho and leads to an estimate of 2.5 kpc. 

Key words: 7-rays: theory - ISM: individual (Tycho) - radiation mechanisms: non- 
thermal 



1 INTRODUCTION 

Energetic cosmic rays (CRs) below the "knee" energy ~ 
3 x 10 15 eV are usually believed to be accelerated in shock 
waves of Galactic SNRs through the diffusive shock accel- 
eration (DSA). Moreover, it has been commonly assumed 
that some 10% of the explosion energy of a supernova is 
transferred to CRs throughout its whole life in order to ex- 
plain the observational CR energy density in our Galax y 
fe.g.. lBaade fc Zwickvlll934lGinzburg fc Svrovatskiilll967f) . 
The theory of DSA of cosmic r ays is well established (see 
e.g.. lMalkov fc O'C Drurvll200ll . for reviews, and references 
therein), but there is no direct evidence for such accelera- 
tion of the CR protons, and the energy conversion efficiency, 
the fraction of the explosion energy converted to CR energy, 
still remain elusive. 

An important clue comes from the hadronic 7-ray emis- 
sion from the neutral pion decay which ensues from the col- 
lision of the accelerated protons with the baryons in prox- 
imate dense matter. It is, however, even difficult to distin- 
guish the hadronic 7-ray emission from the leptonic emis- 
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sion, and thus the environmental effect on the 7-ray emis- 
sion is very noteworthy. So far, 7-rays from some SNRs 
have been interpreted to arise from hadronic processes, such 
as W28, W41, W44, W49B, W5 1C, Cygnus Loo p, IC443, 
CTB 37A, and G349. 7+0.2 (see iLi fc Chenll2012l . and ref- 
erences therein). The 7-rays from the Cas A SNR are also 
explained to have a hadronic contributi on at the GeV en- 
ergy range at least |Arava fc Cuill2010l h The hadronic 7- 
rays from these SNRs are all suggested to be emitted by 
the interaction with dense matter. For the 7-ray emission 
from the Tycho SNR, as a contrast, the h adronic scenario 
invokes a tenuous ambient medium (e.g.. iGiordano et al.,1 
120121 : iTang, Fang fc Zhand[2uTll : iMorlino fc CapriolilbOli T 



Tycho's SNR (G120.1+01.4; 3C 10; SN 1572) is the his- 
torical relic of a Type la supernova as verifie d with the light 
echo from the explosion |Krause et afl feoOSh As one of the 
few accurately known-age SNRs, Tycho's SNR has been well 
studied over the entire spectrum. It shows a shell-like mor- 
phology in radio with enhanced em ission along the north- 
eastern edge (e.g.. |Pickel et al1ll99ll h and has spectral index 
of 0.6 5 and a flux density of 40.5 Jy at 1.4 GHz (|Kothes et al.l 
2006). Strong non-thermal X-ray emission is concentrated 
in the SNR rim, forming a thin filamentary structure (e.g., 
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iHwang et aLll2002l ; iBamba et al.ll2005l ; I Warren et al.ll2005f ). 

Recently, TeV and GeV 7-ray emission was successively de- 
tecte d from this SNR l|Acciari et al, I l201ll ; iGiordano et al.l 
2012) and thus provides a new opportunity for studying 
particle acceleration in SNRs and high energy non-thermal 
radiation. Taking the multi-band spectrum into account, 
it has been shown that the hadronic process, with energy 
conversion efficiency of 10-15% for the curr ent stage, natu- 
rally explains the GeV - TeV 7-ray emission jGiordano et al.l 
120121 ; iTang et alJlioTll ; IMorlino fc Capriolill2012l ). although 
it is argued that a two- zone leptonic model can also work 
(|Atovan fc Der mcr 2012). In the hadronic scenario, the tar- 
get gas is always suggested to be a low-density ambient inter- 
stellar medium (jih ~ 0.2-0.3 cm -3 ), which is quite different 
from the case of the aforesaid SNRs. 

Although it is suggest ed that Tycho is a n aked la SNR 
without any dense clouds jTian fc Leahy|[20TTI ). there seem 
to be signs and evidence in multiband showing the pres- 
ence of dense clouds from the north to the east (see a sum- 
mary in H3.ip . which could be the target for bombardment 
of the accelerated protons and would thus reduce the needed 
converted energy. While, theoretically, the acceleration effi- 
ciency in SNRs is not well determined, very efficient acceler - 
ation process has been sugg ested (e.g. jHelder et al.ll2009h . 
and some self-consistent models predict that the converted 
energy fraction can b e up to ~ 60% - 80% in the whole life- 
time of an SNR fe.g jBerezhko fc Volldli1)971 ; iKang fc Jones! 
2006). However, constraints from some 7-ray bright young 
Galactic SNRs hint that the fraction would not be greater 
than 10% at Tycho's age but could reach 10-20% in subse- 
quent 1-2 kyr, such as the case of RX J1713. 7-3946 (e.g., 
iLee et al.ll2012T ). The similarly young SNR Cas A has only 
transferred a minor fraction (< 2%) of the to t al kinetic en- 
ergy to accelerated particles l|Abdo et al.ll2O10l ; lArava fc Cuil 
|2010T ). Moreover, under the Bohm limit assumption, the con- 
version factor no more than 1.3 % is predicted in kinetic 
models for a SNR at the Tycho's age, evolving in a pre-shock 
medium of density ~ 0.3 cm -3 , althoug h its total conver- 
sion fraction can be more than 60% (e.g.. iBerezhko fc Volkl 
Il997h . Therefore, the fraction of order ~ 1% at Tycho's age 
deserves serious consideration. 

Here we present an alternative hadronic scenario with 
the modest energy conversion fraction of order ~ 1% for the 
present Tycho SNR. As shown below, with such a reduced 
energy fraction, not only the normal number ratio between 
the accelerated electrons and protons is obtained, but also 
the gas density in the region of hadronic process will be 
increased by one order of magnitude, consistent with the 
presence of (at least local) dense medium, e.g, molecular 
clouds (MCs). If this is true, the Tycho's SNR will be the 
only Type la SNR associa ted with MCs amo ng the currently 
known interacting SNRs l| Jiang et al.ll2010l ). 



2 HADRONIC GAMMA-RAY EMISSION 

We assume that the spectrum of the shock-accelerated par- 
ticles (electrons and protons) obeys a power law with an 
energy cutoff: 



dNi/dEi = AiE~ ai x exp(-Ei/Ei,, 



(1) 



where i = e, p, Ei is the particle kinetic energy, cti is the 
spectral index and is fixed as a CjP ~ 2.3 according to the 
afore-mentioned radio index and the charge-independent 
acceleration process, Ei :CUt is the cut-off energy and, for 
protons, can be tak en as -E p , ma x — (37 TeV)(i?d/100/iG) 2 
dParizot et al. 1 120061 ) while -E c ,cut is confined by radiative 
limit, and Ai is normalization factor and can be determined 
from 77-BsN = XT C 2 A c Ec ac x exp(-_E c /_E c , cut )-E c <iE c + 
Im c 2 A P Ep " p x exp( — Ep/E PiCU t)E p dEp combined with ad- 
justable parameter, K ep = A e /A p , the number ratio between 
the accelerated electrons and protons at a given energy. Here 
the energy conversion efficiency n ~ 1% is adopted, and Esn 
is the explosion energy of Tycho. According to the estimate 
from the expansion rate of the X-ray emitting ejecta (|Hughed 
2000), we will adopt E S n = 6 x 10 50 d|, c erg, where cfe.s is 
the distance to the SNR scaled with 2.5 kpc, in view of the 
physical contact of Tycho with interstellar dense clouds (see 
§3). 

For the 7-ray emission, we consider three potential ra- 
diation processes: inverse-Compton (IC) scattering on the 
cosmic microwave background by relativistic electrons, non- 
thermal bremsstrahlung by relativ istic electrons (usin g the 
differential cross sections given in iBaring et al. 19991). and 
7r°-de cay 7-rays resulting from p-p interaction ( Kelner et al.l 
2006, using the parametrization therein). 

The downstream magnetic field strength is estimated 
by X-ray measuremen ts and is listed as a functio n of the 
distance to this SNR i|Cassam-Chenai et ail 120071 . Table 5 
therein). Since in this letter we favor an association of Ty- 
cho with dense clouds at a distance of 2.5 kpc (see §3.3), the 
field strength can be estimated as Bd — 320 /iG by inter- 
polation. This number is also consistent with the estimates 
(around 300 AtG) given b y other methods l|Volk et al.ll2005l : 
IMorlino fc Caprioli|[2012T ). 

From the radio and X-ray synchrotron spectral data and 
the assumed energy conversion efficiency, we get i5 ,cut ~ 
5.6 TeV and K cp « 0.7 x 10 -2 . Actually, the product r/K cp 
is a constant, with the synchrotron flux given (shown in Fig- 
ure[T} ; by contrast, the cases using an efficiency ~ 10% corre- 
spond to a JiTep value < 10 -3 fe.g.. IMorlino fc CapriolHl2012l : 
IGiordano et al.|[2012f l . Recently a lower limit for the electron 
to proton ratio, K ep > 10 -3 , is im plied by the radio obser - 
vations of SNRs in nearby galaxies (|Katz fc W axman 2008). 
Actually, K cp of order ~ 10 -2 is commonly favoured based 
on the measured electron to proton ratio in the CRs under 
the assumption that SNRs are the main source of proton 
and e lectron CRs (see e.g.. lLongairlll99"4l : iKatz fc Waxmanl 
120081) . 

By fitting 7-ray spectrum, which is here considered to 
be due to hadronic emission (see Figure and approx- 
imating the total energy of the relativistic protons to be 
rjEsN, we obtain the product of n and the average density 
(averaged over the entire shock surface) of the target pro- 
tons (with which the energetic particles interact), n t , and 
this product is shown as a n t - n line in Figure [T] We have 
n t ~ 12 cm -3 for r\ ~ 1%. The target protons can be located 
both upstream and downstream, since the accelerated par- 
ticles move back and forth. Downstream, t he average com- 
press ion factor is estimated to be around 3 l|Giordano et al.l 
120121 ). Thus the average preshock target proton density, no, 
is in the range ~ 4-12 cm -3 . This density value is much 
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Figure 1. Dependence of the parameters, K ep and nt, on the 
energy conversion efficiency r\. 



larger than the gas density (nn ~ 0.3cm -3 ), over an order 
of magnitude, estimated from the pre vious optical measure - 
ment for the northeast (NE) filament (Kirshncr ct al. .19871) 
and the Chandra X-ray ob servations (|Cassam-Chenai' et al.l 
l2007l : lKatsuda et al.lfioToh . In £13.21 we will explain such dis- 
crepancy in density for the SNR shock (as in the NE) ex- 
panding in an environment with MCs. 

In the above calculation, the diffusive process of the 
shock accelerated CRs is not considered. Actually, the 
typical diffusion radius of CRs around SNRs is Ran 



2y/D{E)t dit « 0.76 (i dif /440yr) ' 5 pc for 10 GeV protons 
(here the correction factor of slow diffusio n around the SNR , 
Y ~ 0.01, has been adopted, see e.g.. iFuiita et al. 1 120091 ; 
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20121 ). 



iGiuliani et all |2010| ; iGabici et al] I2OI0I ; 

Since the diffusion time scale tan can not be larger than 
the SNR's age 440 yr, the diffusion distance for the protons 
with energy 10 GeV, which are responsible for ~1 GeV 7- 
rays, is smaller than 0.8 pc and significantly smaller than the 
SNR's radius ~ 3 (fe.s pc. Therefore, in the young SNR Ty- 
cho, the hadronic emission does not primarily arise from the 
CR protons which collide nearby molecular gas after a long 
time diffusion, as in the cases of several intera cting SNRs 
that are considerably older (^ 2 kyr, see e.g., iLi fe Chenl 
I2012D . 



3 PHYSICAL CONTACT WITH DENSE 
CLOUDS AND THE DISTANCE 

3.1 The evidence of shock-cloud interaction 

In the previous section, it is shown that a dense ambi- 
ent medium (no ~ 4 - 12 cm -3 on average) is inferred, 
much denser than the preshock gas that was previously 
derived (~ 0.3 cm -3 ). A dense ambient density (3 cm -3 ) 
was also suggested by a recent hydrodynamic simulation 
|f3adenes et al.ll2007l ). Indeed, dense medium (molecular gas 
in particular) has been suggested to be present in the sur- 
roundings of the SNR from the north to the east based on 
multi-band observational evidence. 

The first evidence comes from the reduced expansion 
rate of the S NR in the NE com pared to ot her peripheral 
parts. Radio jRevnoso et al.lll997h and X-ray (Hughes 2000; 
iKatsuda et al. 2010l ) observations show that, along the pe- 



Figure 2. B roadband SEP of T ycho's SNR with the observed 
data in radio jKothes et al . 2006) , X-ray s jTamagawa et alj200sl) 
and 7 -rays (Fermi: lGiordano et al.l2012t VERITAS: I Acciari et aU 
l201lD . 



riphery of the remnant, although the northeastern section 
has the largest radius, its expansion index is smaller than 
that in all the remaining parts. Such evidence suggests that 
the northeastern peripheral gas currently is moving the slow- 
est and being remarkably decelerated and that the gas den- 
sity in the NE is significantly higher than that in other 
parts. It is suggested that the shock front in the NE has 
reached the outskirt s of dense clouds in the recent past 
jRevnoso et al.lll999l ). In addition, the fact that this slowest 
expansion part has the flattest spectral index 0.44 ± 0.02 
ijKatz-Stone et all [2000) mig ht be another indication for 
the interaction with MCs (|Koo et alj l200ll ) . Similar phe- 
nomenon also appe ars in the interactin g regions o f other 
SNR s, such as W28 jDubner et alj|2000l ) and IC443 (|Greenl 
Il986l ). 

Secondly, there seem to be some direct signs of the pres- 
ence of adjacent dense clouds from the CO line observations. 
The 12 CO (J =1-0) observations show that the MCs at 
— 68kms -1 ^ Vlsr ^ — 58kms -1 forms a semi-closed non- 
uniform large shell surrounding the SNR from the north 
to the east a nd the inner boundary essentially traces the 
rim o f Tycho |Lee. Koo fe Tatemats u 2004; Cai. Yang fe Lul 
2009). A high-resolution (16") close-up observation shows 
that the -63.5 km s -1 < VLsr «S -61.5 km s -1 12 CO (J =1- 
0) emitting area appears to be in contact with the remnant 
along its northeastern boundary (|Lee et al.l I2004T ) . Specif- 
ically, there is a big molecular clump in the i ntegration 
map of this velocity component (see Fig. 1 in I Lee et all 
|2004| ) at the azimuth corr esponding to the slowest expan- 
sion noso et all 1 19971) . A similar emitting zone at ~ 
— 62.5kms -1 is also revealed by the submillimeter obser- 
vation in the higher leve l transition line, 12 CO (J =2-1) 
jXu. Wang fe MillerllioTll ). The virial mass of the surround- 
ing MCs derived from the 12 CO (J =1-0) emission is found 
to be larger than their gravitational mass, implying that 
the MCs are being violently disturbed and Tycho' s shock is 
very likely to be responsible for the disturbance (|Cai et al.l 
2009). Such disturbance seems to be consistent with the 
broad 12 CO (J = 2-1) lines pre sent in —68 to — 55kms -1 
interval for a few adjacent MCs |Xu et al.ll201ll ). 

Thirdly, interaction with dense medium is also favoured 
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by the far- to m id-IR observations (|lshihara et alj l20ld ; 
iGomez et alteo^ ). The far-IR emissions (AKARI: 140 and 
160 /J,m; Herschel: 160, 250, 350 and 500 /im), which originate 
from cold dust, are spatially correspondent to the CO emis- 
sion in —68 to — 53kms _1 interval. The mid-IR emissions 
(AKARI: 15, 18, and 24 um; Herschel: 24 and 70 /im), which 
primarily shows warm dust distribution, reveal a shell-like 
emission structure along the remnant boundary except in 
the west and southwest and seem to delineate the interface 
between the remnant and the MCs. They have been pointed 
out to arise from the outer edge of the shock heated cold 
dust and MCs, although the possibility that a small fraction 
of the warm dust was fo rmed in the supernova ejecta was 
not a bsolutely ruled out (jlshihara et alJl20TCt IGomez et al.l 
12012ft . The mid-IR filter bands, such as AKARI L15 (12.6- 
19.4 um) and L18W (13.9-25.6 ^m), may also include the 
contribution from the rotational li nes of hydrogen mole cules 
(like S(l) J = 3-1 at 17.03 /im; llshihara et aI]|201Ch . As 
an alternative, the mid-IR emissions are suggested to come 
from the swept-up interstellar dust through interaction with 
the HI clump in the NE, shown as a significant absorption 
feature in the — 52.7 km s -1 , which is also at the azimuthal 
range of the slowest expan sion and the faint overlapping 
-55kms _1 MC in the NW l| Gomez et alj|2012ft . 

3.2 Discussion about the shock in cloudy medium 

The presence of dense clouds along the boundary, especially 
in the NE, and the large target-particle density derived here 
can be compatible with the low density preshock gas that 
is previously inferred from optical and X-ray observations. 
Here we mention two possible cases for the shock-cloud in- 
teraction. 

The first case is that the blast shock may be flanking a 
cloud (e.g, the slowed-expansion portion in the NE) and the 
shocked gas outside the cloud (but along the line of sight) 
has a density of order 0.1 cm -3 as observed in X-rays and op- 
tical. The solid angle fraction of the whole remnant covered 
by the interacting region may roughly be fa ~ 0.07, where a 
subtend ed angle ~ 60° is adop ted from the slowed-expansion 



lop 

portion l|Revnoso et al.lll997l ). The molecular gas surround- 
ing Tycho may have a density timc as high as ~ 10 2 cm -3 
(T lan fc Leahy||201ir i. Thus the average preshock proton 
density is no = /n(2n MC ) ~ 14(n M c /100 cm- 3 ) cm , 
which is very similar to the value we obtain above (con- 
sidering the uncertainties in the used parameters). 

The alternative case is that the shock may be expand- 
ing into a clumpy medium and the low density corresponds 
to the interclump medium that fills in most of the vol- 
ume. For instance, blast shocks are found to propagate into 
the interclump gas of densit ies of order 0.1 c m" 3 in the 
southeastern edge of Kes 69 iZ hou et alj 2009|) an d in the 
northeastern edge of Kes 78 |Zhou fc ChenTsiOllI ). where 
the shocks are interacting with MCs, as derived from the 
X-ray analyses; in W28, W44, and 3C 391, where blast 
waves are impacting MCs, the preshock density of the in- 
terclump gas is a lso < 1 cm" 3 as infe rred from the IR spec- 
troscopic study (|Reach fc Rholl2000l ). If n ~ 4-12 cm -3 
derived above from the hadronic emission is used and the 
clump filling factor in the dense cloud is denoted as f c (which 
may be of order 0.1 or smaller), then the solid angle frac- 
tion covered by the cloud(s) is fn = io/(2»mc/c) ~ 0.2- 



0.6(riMc/100cm -3 ) -1 (/ c /0.1) -1 . Such a solid angle seems 
to be consistent with the distribution of semi-closed dense 
gas surrounding the SNR as suggested by the CO line and 
IR observations ( i]3.1[l . 

We note that there is not yet report of detection of opti- 
cal emission from the shocked clouds. One of the possibilities 
may be that the shock collision with the clouds is a very re- 
cent event, e.g., within 50 yr (jRevnoso et al.ll"l999t) . which 
is shorter than the r adiative cooling time ~ 300 yr (using 
Eq.(3) in ISgrol 1 19751 . with the velocity of the blast shock 
0.3" yr" 1 at 2.5 kpc or 3500 kms" 1 adopted). Little ther- 
mal X-rays are either detected from the shocked clouds. A 
similar case is seen in RX J1713. 7-3946, anoth er young SNR 
interacting with MCs l|Moriguchi et al.ll2005l ). It is pointed 
out that the substantial suppression of the thermal X-ray 
emission could result from the very low intercloud gas den- 
sity and the stall of the transmitted shock in t he clouds 
( Zirakashvili & Aharonianll2010l : llnoue et alJl2012h . 



3.3 Distance to Tycho 

The SNR-cloud associ ation, in combinatio n with the line- 
of-sight HI absorption l)Tian fc Leahvll201ll ). can be used to 
estimate the distance to the Tycho's SNR, which has long 
been in debate due to its spec ial location in our G alaxy (see 
a brief summary in Figure 6 of lHavato et al . 2010). The LSR 
velocity ~ — 62kms _1 of the semi-closed molecular shell (in- 
cluding the northeastern MCs) allows two alternatives for 
the distance, 4.6 kpc and 2.5 kpc, while the LSR velocity 
~ — 55kms _1 of the faint northwestern MC allows 4.0 kpc 
and 2.5 kpc. The latter estimates are because of the deep 
"dip" at 2.5 kpc, where the spiral shock fr ont is located, 
in th e rotation curve for the longitude 1 20° |Tian fc Leahvl 
l201ll : using the iFoster fc MacWilliamsl ' (|2006h model: also 
see Figure 14 of lSchwarz et allll995l ) and the ^-dependent 
behaviour of the velocity field in the spiral shock region, in 
which the syste matic velocity o f the MCs could be as low as 
< -60kms _1 (|Lee et al.ll2004h . Due to the absence of the 
absorption feature for the HI gas at — 41 to — 46kms~ 1 , Ty - 
cho should be in front of this g iTian fc Leahvl (|201lD 
constraint it to be nearer than 3 kpc, and thus the alter- 
native distances of 4.6 kpc and 4.0 kpc are excluded. In the 
scenario of the association with the ~ — 62kms~ 1 velocity 
compon ent, the HI gas ~ — 60kms~ 1 (lack of absorption 
feature, ITian fc Leahvl l201lf ) is behind the SNR, and the 
—47 to — 53kms _1 HI absorbing gas is in the front, but 
all of them are closely around the spiral shock at the "dip", 
at the distance of 2.5 kpc. In the scenario of the associa- 
tion between the remnant and the ~ — 55kms~ MC and 
the — 52.7 kms -1 HI cloud (see £ ]3.1[1 . the situation of HI ab- 
sorption and SNR location is similar. Tycho is loc ated in the 
high density region just behind the spiral shock l|Lee et al.l 
2004) and is now encountering the inhomogeneous outskirts 
of the ~ -62kms _1 or/and ~ -55kms _1 MCs. 



4 CONCLUSION 

Tycho is one of nearly a dozen Galactic SNRs which are sug- 
gested to radiate hadronic 7-ray emission. It is noted, how- 
ever, that it is the only one in which the hadronic emission 
is proposed to arise from the interaction with low-density 
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(~ 0.3 cm" 3 ) ambient medium. An alternative explanation 
that the 7-rays originate from hadronic process with fewer 
energetic protons but denser target gas is presented here. 
With a conversion efficiency of order 1% for this young his- 
torical SNR, a normal electron-proton ratio (of order 10 -2 ) 
is derived from the radio and X-ray synchrotron spectra and 
an average ambient density that is at least one-order-of- 
magnitude higher is further derived from the hadronic 7-ray 
flux. This result is consistent with the multi-band obser- 
vational evidence of the presence of the surrounding dense 
medium from the north to the east. The interaction scenario, 
combined with the HI absorption data, is used to constrain 
the long controversial problem of the distance to Tycho and 
leads to an estimate of 2.5 kpc. 
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